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THE  ELECTRON  AVALANCHE  IONIZATION  OF  AIR  AND 
SIMPLE  AIR  CHEMISTRY  MODEL 


1.  Introduction 

The  avalanche  ionization  of  air  and  other  gaseous  elements,  caused  by 
a  static  or  an  oscillating  electric  field  (microwaves,  lasers)  is  of 
considerable  interest  in  gaseous  breakdown  studies.  These  studies 
provide  an  understanding  of  the  breakdown  mechanisms  and  facilitate  the 
measurements  of  various  electron  transport  coefficients.  The 
understanding  of  gaseous  breakdown  and  the  ensuing  plasma  chemistry  can 
be  utilized  fruitfully  in  many  areas  of  applications,  such  as,  high 
voltage  switching,  lasers,  plasma  maintenance  and  plasma  heating. 

Avalanche  ionization  occurs  when  a  free  electron  gains  sufficient 
energy,  from  the  externally  imposed  electric  field,  to  generate  an 
ionization,  and  thus  initiate  the  electron  multiplication  process.  The 
initial  free  electron  or  electrons  can  be  introduced  into  the  gas  by 
various  means,  such  as  a  radiation  source,  an  electron  beam,  etc. 

To  describe  avalanche  ionization  in  a  gas,  one  requires  the  solution 
of  the  Boltzmann  equation  with  an  appropriate  set  of  electron-atom  or 
electron-molecule  collision  cross  sections  where  one  obtains,  in  a 
self-consistent  manner,  the  electron  velocity  distribution  and  the 
relevant  transport  coefficients.  However,  the  knowledge  of  several 
quantities  such  as  the  electron  drift  velocity,  the  ionization 
coefficient,  the  electron  temperature  or  the  electron  characteristic 
energy  can  also  be  utilized  to  describe  avalanche  ionization  in  a  speedy 
and  simple  manner.  These  quantities,  in  general,  depend  on  the  electric 
field,  E,  and  the  gas  pressure,  P,  and  have  been  measured**^*^’^  for 
various  gaseous  elements. 
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On  the  other  hand,  the  avalanche  ionization  can  also  be  described  by 
assuming  a  Maxwellian  electron  velocity  distribution;  for  regions  of 
E  -1  -1 

—  >  35  volts  cm  Torr  this  method  provides  identical  results 
compared^  to  the  solution  of  the  Boltzmann  equation. 

In  this  report,  we  discuss  a  simple  description  of  avalanche 

ionization  in  air,  provide  the  ionization  coefficient  over  a  wide  range 
E  E 

of  —  (—  =  50  to  3000)  and  make  comparison  with  previous  results.  The 
simple  approach  consists  of  a  few  equations  which  describe  the  ionization 
and  the  deionization  in  air,  using  appropriate  coefficients  developed  in 
this  report. 

2.  The  Ionization  Coefficient 

The  growth  of  a  current  beam  as  a  function  of  distance,  due  to  the 
avalanche  ionization  in  a  gaseous  medium,  can  be  described  by 

I  »  Iq  Exp  (ox)  (I) 

where  IQ  is  the  initial  current  at  x  =  o  and  o  is  the  ionization 
coefficient  which  is  generally  known^  as  Townsends'  first  ionization 
coefficient.  This  coefficient,  especially  in  the  form  of  has  been 

E 

measured  as  a  function  of  —  for  various  gaseous  elements  results  are 
quoted  in  Ref.  5,  but  in  that  reference  reliable  data  for  air  and  its 
constituents  are  limited  to  the  region  of  E/P  below  100  volt  -cm~* 
-Torr-1.  These  units  of  —  are  used  throughout  this  report  and 

O 

henceforth  implied  throughout  this  report.  The  work  of  Rao  and  Raju  , 

E 

however,  extends  the  range  to  »  1000  for  dry  air  with  no  contribution 

from  contaminants.  Their  data  for  air  is  given  in  Table  1  for 
E  E 

—  *  54  to  *  1000  and  is  also  shown  in  Figure  la.  The  data  of  Table  1 
can  be  expressed  reasonably  well  (within  few  percent)  by  the  following 

O 

expressions. 
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f  -  8.34  Exp  (-  |^i),  120  >|  >  54  (2) 

and 

1000  >|  >  120 

(3) 


Table  1 


The  Mr 

Ionization  Coefficient 

(p  ) As  a  Function  of  (E/P) 

E/P 

0  /P 

E/P  o/P 

54 

5.24  (-2) 

330 

4.95 

60 

8.70  (-2) 

360 

5.47 

66 

1.32  (-1) 

390 

5.96 

72 

1.86  (-1) 

420 

6.41 

80 

2.72  (-1) 

450 

6.82 

90 

3.98  (-1) 

480 

7.21 

100 

5.40  (-1) 

520 

7.68 

120 

8.52  (-1) 

540 

7.9 

130 

9.93  (-1) 

580 

8.69 

140 

1.17 

620 

9.11 

150 

1.36 

660 

9.49 

160 

1.55 

700 

9.83 

170 

1.73 

740 

10.15 

180 

1.92 

780 

10.45 

196 

2.20 

820 

10.73 

210 

2.49 

860 

10.98 

230 

2.94 

900 

11.22 

250 

3.37 

940 

11.44 

270 

3.79 

980 

11.65 

300 

4.39 

Q 

The  ionization  coefficient  for  N2  as  measured  by  Folkard  and  Haydon 
for  E/P  from  60  to  1000,  and  that  for  O2  as  measured  by  Price,  et  al^for 

g 

from  30  to  60  are  shown  In  Figure  lb.  The  extension  of  the  ionization 
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coefficient  for  air  and  N2  down  to  E/P  ■  40  is  based  on  data  from  Ref. 
5.  However,  Equation  (2),  if  used  below  its  lower  limit  will 
underestimate  ^  by  ~30%. 

3.  The  Drift  Velocity 

The  electron  drift  velocity  in  N2  and  O2  has  been  measured  by 

11  p  r 

Schlumbohnr 1  over  a  wide  range  of  ^  .  For  N2  the  range  of  —  is  from 
100  to  3000,  while  for  O2,  the  range  is  from  100  to  8000.  Using  these 
data**,  the  drift  velocity  for  N2  and  O2  can  be  expressed**,  reasonably 
well,  by  the  following  relations 


v 


d 


120  <  |-  <  3000  (4) 


v 


d 


3.75  x  106 


100  <  |-  <  8000  (5) 


respectively,  where  the  units  of  v^  are  in  cm  see-*.  The  drift  velocity 
in  air  can  be  obtained  using  the  relation^1  for  the  electron  mobility,  p  , 
in  a  gaseous  mixture  where 


(6) 


Here,  fj  and  f2  are  the  fractions  of  N2  and  O2  in  air  and  p^  and  p^  are 
the  electron  mobility  in  N2  and  O2,  respectively.  From  Eqs.(4-6),  we 
obtain  the  following  expression  for  the  electron  drift  velocity  in  air 


v 


d 


-  3.38  x  106 


5 


100  <  |-  <  3000  (7) 


'  n 


r 


- -7- 


For  the  region  of  -^  <  100  the  electron  drift  veocity  in  N£  and  0£ 

1  2 

are  shown  in  Fig.  2  based  on  the  data  of  Wagner  and  Raether  ,  and  Jager 
and  Otto^.  The  measurements  of  References  12  and  13  with  Equation  (6) 
predict  the  drift  velocity  in  air  which  is  in  good  agreement  with  the 
drift  velocity  measurements  of  Ryzko*\  shown  in  Fig.  2,  and  is  expressed 


vd  *  6.0  x  106  +  2.5  x  10s  (E/P)  (8) 

1  C 

Felsenthal  and  Proud  have  proposed  the  following  expression  for 

the  electron  drift  velocity  in  air 

vd  *  7  x  106  +  2  x  105  (|  )  (9) 

£ 

based  on  measurements  over  a  limited  range  of  E/P  (—  *  0.2-20)  .  They, 

£ 

however,  extended  Equation  (9)  to  higher  values  of  —  which  clearly 

overestimates  the  drift  velocity.  The  dependence  of  the  drift  velocity 

F  E  E  ^2  £ 

on  is  linear  for  <  100,  and  varies  as  (—  )  for  —  >  100  as  shown 

£ 

by  Equations  (7)  and  (8).  For  example  at  -p  *  1000  Equation  (9)  predicts 

ft  O 

a  drift  velocity  of  2.07  xi0°  cm/sec  compared  to  1.06  x  10  cm/sec  based 
on  Eq.  (7)  which  is  obtained  from  the  measurements^  of  the  drift 
velocities  in  N2  and  O2. 

4.  The  Ionization  Frequency  in  Air 

The  ionization  frequency  is  related  to  the  ionization  coefficient 


Vi  "  °  Vd 


or  alternately 


'  ^  * 


7 


(11) 


p 


a 

P  Vd 


Vi  E 

Thus,  — can  be  expressed  as  a  function  of  —  through  the  dependence 

of  both  —  and  ^  on  —  .  Using  Equations  (7)  and  (8)  for  the  drift 

velocity  and  Equations  (2)  and  (3)  for  we  obtain  the  following 

expressions  for  the  ionization  frequency  in  air. 


*  1.0  xlO7 [5  +  0.21  (E/P)]  Exp  (-  -^7^)  .40  s<|  s<100  (12) 


and 


~J  «  4.94  x  107  (  |  )1;2  Exp  (-  1^1  ) 


,  120  x<  |  x<1000 


(13) 


i  E 

Figure  3  shows  —  as  a  function  of  —  based  on  Equations  (12)  and  (13). 

However,  the  ionization  frequency  can  also  be  obtained  from  the 
basic  atomic  processes  if  the  electron  velocity  distribution  is  known. 

For  a  Maxwellian  electron  velocity  distribution,  the  ionization  rate 
coefficients  for  O2  and  N2  have  been  calculated  *7  using  measured 
ionization  cross  sections.  These  rate  coefficients  can  be  used  to  obtain 
the  corresponding  ionization  rate  coefficient  in  air.  The  numerical 
values  for  the  rate  coefficients  in  ??2  and  C>2,  shown  in  Figure  4,  can  be 
expressed**  as 


.9 

<<JV>„  -  (-4.98  +  9.67  I  +  0.49  Ta2)  Exp  (-  15‘ -  ) 


S2 


e 
-9 

<ov>  -  (-2.02  +  3.39  T  +  0.64  T  2)  Exp  (-  y-‘°-) 

°2  fre  e  e  Te 


(14) 


(15) 


’  U 


These  expressions  are  valid  for  electron  temperatures  from  0.5  to  15  eV 
and  deviate  by  less  than  10%  from  the  manerical  values*7.  Using 
Equations  (14)  and  (15)  the  ionization  frequency  in  air  can  be  expressed, 


~  -  12_  (14.093  +  27.366  T  +  1.386  T  2j  Exp  (-  ) 

ff.  e 

+  —  (1.430  +  2.470  T  +  0.456  T  2]  Exp  (-  y- °--  )  (16) 

/"ST  e  e  T 

>Pe 


where  Te  is  in  units  of  eV.  The  ionization  frequency  in  air  (Eq.  16)  is 

shown  in  Figure  5  as  a  function  of  the  electron  temperature. 

1 8 

Lupan  has  developed  an  expression  similar  to  Equation  (16)  based 

on  a  Maxwellian  electron  velocity  distribution  and  assuming  that  the  main 

1  8 

ionization  is  from  02*  His  expression  reads 

~  *  2.5  x  107  (5.6  (E)1/2  -f  0.6  (E)3/2  ]  Exp  (  -  18 5  )  (17) 

P  E 

where  E  ■  1.5  Tg  .  Equation  (17)  predicts  consistently  higher  values  for 
the  ionization  frequency  compared  to  Eq.  (16)  .  The  discrepancy  is  a 
factor  of  2  at  Te  *  2eV  and  14%  at  Te  =  lOeV,  and  is  clearly  due  to  the 
assunption  by  Lupan  that  ionization  in  air  is  mainly  from  02»  A  glance 
at  Figure  4  shows  that  the  ionization  rate  coefficient  in  O2  is  always 
higher  at  lower  Tg  compared  to  that  for  air.  At  higher  Tg  the  two 
converge  to  the  same  value,  indicating  that  at  higher  temperatures  the 
ionization  rate  coefficient  for  N2,  O2  and  air  are  practically  equal. 
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5 .  The  Electron  Temperature 

The  electron  temperature,  in  an  electric  discharge,  can  be  related 
E 

to  —  .  In  general,  however,  one  defines  an  electron  characteristic 
i  2 

energy  ’  ,  ,  which  for  a  Maxwellian  electron  distribution  is 

equivalent  to  the  average  temperature,  i.e.,  «  Tg  . 

In  this  work  we  would  like  to  obtain  the  electron  temperature  in  air 
as  a  function  of  E/P  by  utilizing  the  ionization  frequency  which  was 
obtained,  as  a  function  of  E/P  and  as  a  function  of  Te  (see  Figures  3  and 
5,  respectively).  From  the  results  of  Figures  3  and  5,  one  obtains  a  one 
to  one  correspondence  between  Tg  and  (E/P).  This  correspondence  is 
expressed  as 


where  Tg  is  in  eV.  Eq.  (18)  gives  a  good  correspondence  to  better  than 
15%  for  E/P  from  60  to  500.  At  p—  “  1000  ,  however.  Equation  (18) 
overestimates  the  temperature  by  22%.  Equation  (18)  is  shown  in  Figure  6 
and  is  compared  with  various  results  There  is  generally  a 

good  agreement  between  Eq .  (18)  and  the  results  of  Reference  2  and  6 
where  the  data  is  for  N2  and  it  is  expected  that  the  electron  temperature 
is  somewhat  higher  in  N2  than  in  air.  Lupan' s  result  ,  valid  for 
Te  <15eV,  is  always  higher  than  our  calculation  as  given  by  Eq.  (18). 
Therefore,  the  extension  of  Lupan' sl  result  to  -p  >  100  will  always 
overestimate  the  electron  temperature. 
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6.  Avalanche  Ionization  in  Air 

To  describe  Che  ionizaclon  in  air,  a  sec  of  race  equacions  which 
gives  Che  cime  hisCories  of  various  ionic  and  air  species  is  required. 
These  race  equacions  muse  be  solved  simulcaneously  wich  an  equacion 
describing  Che  elecCron  energy  or  Che  elecCron  cemperaCure.  However,  for 
simplicicy  fewer  race  equacions  can  be  udlized  Co  generaCe  a  reasonable 
descripCion  for  Che  avalanche  generaced  ionizaCion.  In  Chis  secCion,  we 
provide  such  a  descripCion  wich  Che  appropriaCe  parameCers. 

We  consider  Che  case  of  dry  air,  however,  Che  ionizaCion  of  air  wich 
impuricies,  e.g.,  HjO  will  be  discussed  in  a  fuCure  reporC.  The 
equacions  describing  Che  avalanche  ionizaCion  can  be  wriCCen  as 


dT 
e  dc 


e 


■<V7  V 

s  L  e 


-  He02R  (02) 


S  +  Y  P 


2‘\ 


mv 


m 


+  NeHv 


-  Ne  N2  R  (N2) 
(19) 


dN 
_ e 

dc 


S  +  N 

e 


ef 


o  N  0«M  -a.  N  0. 
a  e  2  da  e  2 


(20) 


dN 

dC 


a  N 
a  e 


02  M 


+  a 


da 


N  0, 
e  i 


a  N_(N  4  N  ) 
m  e 


(21) 


dE  8  vE. (T  -T  ) 

“  NeE1  (1-Exp(-E1/Ty);  l  vXyU-Exp(  -  -  )]  (22) 

v»l  v  e 
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Here  S  is  the  volume  production  rate  (cm  sec  )  of  electrons,  due  to  any 

other  source,  £  is  the  average  energy  of  these  electrons,  y  is  a 
s 

conversion  factor,  Ng  is  the  plasma  electron  density,  is  the  momentum 

transfer  collision  frequency  (electron  neutral  plus  electron  ion),  R(^) 

is  the  rate  of  energy  loss  by  the  plasma  electrons  in  N2  and  R(02)  is  the 

corresponding  rate  in  O2.  <*e£  is  the  effective  recombination  rate 

coefficient,  o  is  the  coefficient  for  electron  attachment  and  a,  is 
’a  da 

the  dissociative  attachment  rate  coefficient.  a  m  is  the  coefficient 
for  the  mutual  neutralization.  N  is  the  sum  of  the  negative  ions,  and 
that  the  charge  neutrality  condition  is  assuned.  Ey  is  the  vibrational 
energy^  0£  jj  £s  the  energy  for  the  N2  vibrational  spacing  (assumed 

a  harmonic  oscillator) ,  Tv  is  the  vibrational  temperature,  Xy  is  the 
excitation  rate^  coefficient  for  the  v  —  vibrational  level  and  Hy  is 
the  vibrational  heating  rate  of  the  electrons  and  corresponds  to  the  sum 
of  the  negative  terms  in  Eq.  (22). 

The  rates  of  the  energy  loss  by  electrons  in  N2  and  C>2  as  a  function 
of  the  electron  temperature  are  shown  in  Figures  7  and  8,  respectively. 
These  are  obtained  by  using  the  relevant  rate  coefficients^  and  the 
appropriate  threshold  energies^.  In  these  figures  the  energy  loss  to 
various  processes  are  also  indicated.  These  rates  can  be  expressed  as 

R(N2)  -  3.87  x  lCf8  Tg  Exp(-  |^-)  ,0.2$T^1.0  (23) 

e 

R(N2)  -  7.5  x  icf9  tl-0.036Tel  +  A  ,1$T,<20  (24) 
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where 


A  -  2.13  x  10"7  pT  [l  -  0.05  (T  -4)]  Exp  (-  )  (25) 

’  e  e  r 

e 


and 


R(02)  -  1.0  x  10 ~7fT  [1  +  0.05  (Te-4)]  Exp  (-  )  (25) 

e 


In  addition  to  Equations  (19-22),  one  oust  solve  for  the  time  variance  of 
the  electric  field  which  enters  into  Equation  (19).  This  can  be 
accomplished  by  solving  an  appropriate  electric  circuit  equation. 

The  effective  recombination  coefficient,  ,  can  be  obtained, 
during  the  discharge,  by  appropriate  steady  state  relations  for  the 
density  of  various  ions.  These  relations  must  account  for  all  source  and 
depletion  terms  for  each  ion.  For  the  case  of  the  dry  air  with  no  air 
ion  clusters,  the  effective  recombination  coefficient  can  be  expressed  as 


No 
e  ef 


D!^(N2+) 


DlNe  +  C2  °2  +  V 


d2q(02+) 

D-N  +  o  N_ 
z  e  m 


d2c2o2  Q(N2+) 

(D.N  +  a  N~)  (D. N  +  a  N~+  C,0») 
Ze  m  le  m  z  z 


DjCj  02  Q  (N  ) 


(D.N  +  a  N  )  (C,0  +  o  N  ) 
z  e  m  l  z  o 


(27) 


where  Q(i+)  is  the  volune  production  rate  of  the  1^—  ion  due  to  the 
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ENERGY  LOSS  RATE  COEFFICIENT  IN  N2  (eV*cm3-sec“ 


ENERGY 


Fig.  8  -  The  electron  energy  loss  rste  in  Oj  as  a  function  of 


source  Cera  and  Che  avalanche  Ionization.  Dj  and  D£  are  Che  dissociative 
recombination  rate  coefficients  of  Nj+  and  0^+  ,  respectively;  and  Cj  and 
C2  are  the  charge  exchange  rate  coefficients  of  N+  and  Nj+  with  O2, 
respectively.  These  coefficients  can  be  found  in  Reference  21.  However, 
if  the  steady  state  relations  overestimate  the  positive  ion  densities 
especially  at  early  times,  one  can  use  the  following  relation 

“ef  Ne  *  {K1D1+  K2D2  +  **J  *Ne  +  (28) 

where  Kj  and  K2  Indicate  the  fractions  of  the  appropriate  ions  in  the 
discharge. 
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